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Q ' 1. Introduction 

' The geodesic flow on a Riemannian manifold represents the extremals of the least action 

principle, namely it is determined by the motion of a certain physical system in the manifold. 
It is known that the geodesic equations are second order non-linear differential equations and 
they usually appear in the form of Euler-Lagrange equations of motion. Magnetic curves 
generalize geodesies. In physics, such a curve represents a trajectory of a charged particle 
moving on the manifold under the action of the magnetic field. 

Let (M, g) be an n-dimensional Riemannian manifold. A magnetic field is a closed 2-form F 
\ on M and the Lorentz force of a magnetic field F on (M, g) is an (1,1) tensor field 3> given 

cn ■ b y 

9 mX),Y)=F(X,Y), VI,F6x(M). (1) 

The magnetic trajectories of F are curves 7 on M that satisfy the Lorentz equation (some- 
times called the Newton equation) 

Vy 7 ' = $(</). (2) 
, Lorentz equation generalizes the equation satisfied by the geodesies of M, namely 

S : v y y = 0. 

Therefore, from the point of view of the dynamical systems, a geodesic corresponds to a 
trajectory of a particle without an action of a magnetic field, while a magnetic trajectory is 
a flowline of the dynamical system, associated with the magnetic field. 

Since the Lorentz force is skew symmetric we have 

^ ff ( 7 / ,7') = 2 ff (V r7 ',7 , ) = 0, 

so the magnetic curves (trajectories) have constant speed v(t) = || 7 '|| = vq. When the 
magnetic curve j(t) is arc length parametrized (vq = 1), it is called a normal magnetic 
curve. 

Recall that a vector field V on M is Killing if and only if it satisfies the Killing equation: 

g(V Y V,Z)+g(V z V,Y)=0 (3) 
for every vector fields Y, Z on M, where V is the Levi Civita connection on M. 
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A typical example of uniform magnetic fields is obtained by multiplying the volume form 
on a Ricmannian surface by a scalar s (usually called strength). When the surface is of 
constant Gaussian curvature K, trajectories of such magnetic fields are well known. More 
precisely, on the sphere E> 2 (K), K > 0, trajectories are small (Euclidean) circles of radius 
(s 2 + K)^ 1 / 2 , on the Euclidean plane they are circles and the period of motion equals to — , 
while, on a hyperbolic plane H 2 (— K), K > 0, trajectories can be either closed curves (when 
\s\ > Vif), or open curves. Moreover, when |s| = \f~K normal trajectories are horocycles 
(see e.g. [IIJEE2]). 

This problem was extended also for different ambient spaces. For example, if the ambient 
is a complex space form, Kahler magnetic fields are studied (see [2]), in particular, ex- 
plicit trajectories for Kahler magnetic fields are found in the complex projective space CP™ 
PQ. Kahler magnetic fields appear in theoretical and mathematical physics, varying from 
quantum field theory and string theory to general relativity. 

If the ambient is a contact manifold, the fundamental 2-form defines the so-called contact 
magnetic field. Interesting results are obtained when the manifold is Sasakian, namely the 
angle between the velocity of a normal magnetic curve and the Reeb vector field is constant 
(see 0). Moreover, explicit description for normal flowlines of the contact magnetic field on 
a 3-dimensional Sasakian manifold is known 

In the case of a 3-dimensional Riemannian manifold {M,g), 2-forms and vector fields may 
be identified via the Hodge star operator ★ and the volume form dv g of the manifold. Thus, 
magnetic fields mean divergence free vector fields (see e.g. [10]). In particular, Killing 
vector fields define an important class of magnetic fields, called Killing magnetic fields. It 
is known that geodesies can be defined as extremal curves for the action energy functional. 
A variational approach to describe Killing magnetic flows in spaces of constant curvature is 
given in [7]. 

Note that, one can define on M the cross product of two vector fields X, Y £ x(M) as follows 

g(XxY,Z) = dv g (X,Y,Z), VZ e X (M). 

If V is a Killing vector field on M, let Fy — ivdv g be the corresponding Killing magnetic 
field. By t we denote the inner product. Then, the Lorentz force of Fy is (see [TU] ) 

$(X) = V x X. 

Consequently, the Lorentz force equation ([2]) can be written as 

V 7 / 7 ' = V X y. 

In what follows we consider the 3-dimensional Euclidian space E 3 , endowed with the usual 
scalar product ( , ) . 

The fundamental solutions of ^ are {d x , d y , d z , —yd x + xd v , —zd y + yd z , zd x — xd z } and 
they give a basis of Killing vector fields on E 3 . Here x,y, z denote the global coordinates on 
E 3 and R 3 = spanfc^, d y ,d z } is regarded as a vector space. 

The easiest example is to consider the Killing vector field £o = d z . (Similar discussions 
can be made for d x and d y , respectively.) Its trajectories are helices with axis d z , namely 
1 1 — ^ (xq + a cost, yo + asint, zo + bt), where (xq, yo,zo) g M 3 and a, b € K. An interesting fact 
is that Lancret curves (i.e. general helices) in E 3 are characterized by the following property 
(in our framework): they are magnetic trajectories associated with magnetic fields parallel 
to their axis. A similar result, relating Killing magnetic fields and Lancret curves is provided 
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on the 3-sphere (see e.g. [7]). Theorems of Lancret for general helices in 3-dimensional real 
space forms are presented in [3]. 

In this paper we consider the following magnetic field Fy = —(xdx + ydy) A dz in E 3 , 
determined by the Killing vector field V = —yd x +xd y . The other two rotational vector fields 
—zd y + yd z and zd x — xd z give rise to analogue classifications for corresponding magnetic 
trajectories. The aim of this note is to find all magnetic curves corresponding to Fy. The 
main result we obtain is the following: 

Theorem. The magnetic trajectories of the Killing magnetic field Fy are: (a) planar curves 
situated in a vertical strip; (b) circular helices and (c) curve parametrized by 

t 

x(t) = p(t) cos c/>(t), y{t) = p(t)smct>(t), z(t) = ~\ J P 2 (CR 
where p and <f> satisfy 

("dP) + P ( p2(<) ) = °' ^ WW = constant 

and P is a polynomial of degree 3. 

We are able to obtain explicit solutions in case (c) and we represent some examples by using 
numerical approximations for some integrals. 

Recall, for later use, some basic facts on normal elliptic integral of the first kind (see for 
example [5]): 

y ip 

U = sn _1 (y, k) = F(ip, k), 



^(l-t^(l-kH^ J Vl~fc 2 sin 2 tf 
u o 

where y — sirup and ip = am u. The angle ip is called Jacobi amplitude and the function 
sn in known as Jacobi elliptic sine. The number k is called modulus and for applications to 
engineering and physics it belongs to (0, 1). 



2. Rotational magnetic trajectories in E 3 

Let us consider the Killing vector field V = —yd x + xd y on E3 \ Oz, which defines the 
magnetic field Fy = —(xdx + ydy) A dz. The Lorentz force <&y acts on the vector space R 3 
as follows: 

®vd x = -xd z , <&yd y = -yd z , § v d z = xd x + yd y . 
For the Euclidian space E 3 the Lorentz force equation becomes 

7" = V X 7' (4) 

where the curve 7 : J = [0,Z] — > E 3 , j(t) — (x(t),y(t), z(t)) is parametrized by arc length, 
namely 

x'(tf +y'(tf + z\tf = 1, Vie/ (5) 

and at the moment t = it passes through the point (xo, yo, z$), with the velocity (uq, vq,wq), 
such that 

U Q+ V 0+ W = 1- 
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Proof of the Theorem. Our aim is to determine the magnetic curves of Fy- The equation 
(jlj yields the following ordinary differential equations system 

x" — xz' 

y" = yz' (6) 

y = -{xx' + yy'). 

In order to solve it, note that from the first two equations we get a prime integral 

x'y - y'x = u y - x v (7) 
while from the third equation we obtain 

^' = -\{x 2 + y 2 ) + \{xl + yl) + w . (8) 

Notice that z' cannot vanish identically (on a subinterval of I). Indeed, if z' = then 
x' = Uo, y' = v and z — z with u 2 , + v 2 , = 1. Hence, x(t) = Xq + u t, y(t) — y Q + v t 
and combining with © we get a contradiction. It follows that one cannot have horizontal 
magnetic curves corresponding to V. 

In the sequel it is more convenient to consider cylindrical coordinates {p, 4>, z} on E 3 \ Oz. 
Thus, for our curve we have 

x = p(t) cos 4>(t) 
V = P{t) sin 
z = z(t) 

where p 2 (t) = x 2 {t) + y 2 (t), p(t) > 0. 

Case I. First we study the general case, when z' is not constant (equivalently p is not 
constant). The relations ([7]) and ([8]) lead to 

p 2 (t)0'(t)= Po (9) 

z'(t)=qo-\ P 2 {t) (10) 
where we put p Q = x Q v - u y and q Q = \ (x\ + yl) + w - 
The arc length parametrization condition (J5J) , together with (|10p . becomes 

p' 2 (t) +p 2 (t)4>' 2 (t) +q 2 Q - q p 2 (t) + \p\t) = 1. (11) 

Multiplying (|TTJ) by 4 j o 2 (t), using ((9]) and denoting p 2 (t) by f(t) > 0, for all t e I, one gets 

!' 2 + f- 4<Zo/ 2 + 4(q 2 - 1)/ + 4p 2 = 0. (12) 

We start to study the above differential equation for some particular values of the constants 
po and go- 
lf po = 0, i.e. xqVq = youo it follows that the angle <fi is constant, <f> — 0o, so the magnetic 
trajectory is a planar curve, with 

x(t) = p(t) cos(j> , y(t) = p(t) sin (f> . 

More precisely, the curve lies in the plane (suk/)q)x — (cos</>o)y — 0. The initial conditions 
expressed in cylindrical coordinates, may be written as 

.t = po cos 4> , y = Po sin </> 
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and the condition xqvq = youo becomes uo = Cocos</)o, vq — (bsin<^o, for a certain Co G R. 
It follows that Co + (<7o — | Po) = 1 

-1 + ^ Po < Qo < 1 + ^ Pi 
Since po ^ it follows that go > — 1- 

Let us solve the equation (|12p. for three particular situations arising from the initial condi- 
tions: 

(i) If pq = and go = 0, then the equation (fT2l) takes the form 
f 2 (t) + f(t)(f(t)-2)(f(t) + 2)=0 

and it has solution if and only if f(t) < 2, i.e. p(t) 6 (0, \/%, so the magnetic curve 
7 lies inside a cylinder. In fact, being a planar curve, 7 stays in a vertical strip 
centered in Oz and of width 2v / 2. 

We have /'(*) = ~^*J f(t) (4 — fit) 2 ) and we consider only the plus sign (the other 
situation may be treated in similar way). Supposing po ^ ^/2, we have that / and 

the integral 1(f)— / — — are strictly increasing functions. Thus, the 

hi VC(4-C 2 ) 

equation = t has a unique solution in the interval (pg,2), namely / = J{t), 
where J is the inverse function of X. Consequently, p(t) = \J 3{t). In fact J may 
be expressed in terms of the elliptic functions. More precisely, 

2sn 2 (t + t ,^) 
J{t) = 2-sn^ + to,-^) 

where to is determined by sn(£o, ^) = J^° 2 ■ 
Summarizing, the magnetic curve is given by 



*(*) = V^Wcos^o, y(t) = V^Wsin^o, «(t) = -g j o ^(CR- 

In order to draw a picture of our curve, one can use Matlab to compute the 
parametrization. The idea is to calculate the integrals numerically, as Riemann 
sums. See Appendix. 

(ii) If po = 0, g = 1, then the equation jT2]) becomes f' 2 (t) + f 2 (t)(f{t) - 4) = 0, from 
which we have that f(t) < 4, equivalently p(t) < 2, so the magnetic curve 7 stays 
inside a cylinder of radius 2. In fact, being planar, the curve lies in a vertical strip 
centered on z-axis. The equation can be written in the form 

df 



Taking the plus sign, one gets the solution 

4 



± dt. 



/(*) = u2u — , *e(o,t ) 

cosh (t — to) 



where t = — i In 2 p ° . Hence 

2 2+^/A-pl 



P ^ cosh(t - t ) 
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and the magnetic curve is parametrized by 

. . 2 cos 4>o , s 2 sin (j>o „/■•■/ s \ 

x ( t ) = m — n> 2/W = m — TT' z \ t ) = zo + t-2 (taring - to) + tanhi ) . 

cosh(i — to) cosh(i — t ) 

We draw a picture of this (planar) curve. 



Case p Q =0 and q Q =1 




Figure 1. t G 



Let us finalize the examination of the equation (jf 2[) for po = 0. The polynomial 
P(f) = /(7 2 -4<?o/+4(^-l)) has three real solutions, namely fx = 2(q -l), f 2 = 2(q + l) 
and fz = 0. If / is a solution for (fT2")l . then P(/) should be negative. Recall that qo > — 1. 



We have 



(a) If q G (-1, 1), then fx < < / 2 . It follows that p(t) € (0, ^/2(q + 1)) and the 
discussion is similar as in case qo = 0. More precisely we have 

2(l- gn 2 )sn 2 (< + to, V / i ^ 1 ) 

= 



(<7q + 1) six 2 (t + t ,^ 



where to is defined by sn(to, 



22+1 \ _ 2 Po_ 



2 I V * + 1 Vpo-2(«o-1)' 

(b) If (jo > 1, then < f\ < f 2 . It follows that f(t) G (/i,/^)- Thus, the curve 7 lies 
between two cylinders since G ^\/2(<7o — 1), \/2(<Zo + 1)) • As before, the curve 
is situated in a union of two vertical strips. Again, the discussion is similar as in 
case qo — 0. In terms of elliptic functions, we may write 



P(t) 2 " *" 



22 ±l_ sn 2 (v / 2a ±l t + t0) ^_2_) 

where to is defined by sn(t , yfg^) = 

In order to visualize an example, consider the following initial conditions: xo = 2, 
yo = 0, zo = and uo = 0, «o = 0, wo = 1 (this yields po = and go = 3)- 
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We will use again Matlab to compute the integrals (numerically) and to draw the 
picture. 



Magnetic curve in E 3 




Figure 2. p = 0, q = 3 

Return to (p~2]) for po ^ and notice that the equation 

P(f) = f - 4?o/ 2 + % 2 - 1)/ + Ap 2 = (13) 

has the discriminant 

A = -16 [27^ + 8p 2 q (q 2 - 9) - 16( 9o 2 - l) 2 ] 
and the following situations appear: 

• the equation (fl"3|) has three distinct solutions iff A > 0. 

• the polynomial P has multiple roots iff A = 0. 

• the polynomial P has one real root and two complex conjugate roots iff A < 0. 

A detailed analysis of the above situations, lead us to conclude, after taking into account 
classical Viete's formulas, that the equation (fT2j) has solutions if and only if A > 0. 

Indeed, if A < 0, let A E C \ K and A be the complex solutions of (fT3"]) , and B its real 
solution. Then, the ODE (fT2|) can be rewritten as 

f'(tf + (f(t) 2 - 2 Re(A) t + \A\ 2 ) (f(t) -B)=Q, 

where Re(A) denotes the real part of the complex number A. From the third Viete's formula 
we conclude that B should be negative, and consequently, the previous equality cannot occur. 

On the other hand, if A — 0, analyzing the coefficients one cannot have a triple root (since 
Wq 2 ^ 12(go - 1). Hence, let A e R be the double root, and let B G M be the third one. 
With a similar argument as above, B is negative and the ODE (fT2")l becomes 

f(t) 2 + (f(t) A) 2 (f(t) -B)=Q. 
Again, this equality cannot hold. 

It follows that A should be (strictly) positive. Let A, B, C £ R be the three distinct 
solutions of (|13| . The third Viete's formula yields ABC = — \p\ < 0, and hence 
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a) either A, B, C are all negative, 

b) or two of them, A and B, are positive and the third one, C, is negative. 

In case a) the ODE CE2]): f'(t) 2 + (f(t) - A){f{t) - B){f{t) - C) = has no solution. This 
happens if and only if qo < — 1 and po ^ (for the proof use the second and the third Viete's 
formulas) together with A > (for example if qq = — 3 and po = 1). 

In case b), equivalently to A > 0, qg > — 1 an d Po 7^ 0, the equation (JTSJ) has a solu- 
tion in the interval defined by the positive solutions A < B of (|13p . Since the function 

I(f)= / — , is strictly increasing, 1(f) = t has a unique solution 

Ja ^(C-A)(B-C)(t-C) 

r* d( 

f, denoted by J{t). Thus we have p = \J J{t), and <f>{t) — </>o + Po / , , . In this case, 

Jo o[Q) 

the magnetic curve 7 is given by 

x(t) = y/j(ficoa(4>o+PoJ jffi), V® = V^V) wo. +Po J jkr), 

1 /■' 

z(t) = z + q t-- / J(()d(. 

1 Jo 

We may express p in terms of elliptic functions, namely 

Pit? = 



, _ Ak 2 - Csn 2 (rt + t Q , ±) 



fc 2 -sn 2 (rt + i ,£) 



where k 2 = |=§, r - and an(t 0> £) = fc^J 



-c 

In the Appendix we will draw a picture (using the same technique in Matlab as before) 
corresponding to the following data: p — ^/ 2v ^~ 3 ; g Q — 3 ~>^ , for which we have A = 1, 
B = 2 and C = 3 — 2-y/6. 

The situation /(/) = — t can be treated in similar way. 

Finally, notice that for q n = —1 we get A = — 16pq(27j>q + 64) and this case was discussed 
above. 

Case II. Now, let us study the remaining case when z'(t) = wq 0. We immediately have 
that 

z(t) = z Q + tw , 

and from ijHJ) we obtain 

x 2 + y 2 = x 2 + y 2 . (14) 
This means that the magnetic trajectory 7 lies on the circular cylinder of radius 
Po = x/xl + Vo- 

Two subcases must be discussed: Wq < and Wq > 0. 

II. 1: In the case when wq < the magnetic curve is given by 
x(t) = x cos(^-w t) + -j= sin(V-wo t) 

y(t) = y cos{y/^t) + ^^sm(y/=WG~t) (15) 

z(t) = Zq + tw . 

This curve is a helix around the above cylinder. 
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At this point, we have to find which are the initial conditions leading this situation. 
To do this, using (|14[) and (fT5|) . we should have the following relations 

2 . , 

wo = — , , , wo = £pov-wosin</>o, v = -epoV-wocos^o 

Pi + 

where p and </>o have the usual meaning and e = ±1. 
II. 2: If wq > 0, the ODE system ^ has the following solution 

x(t) = x Q cosh(^/wo" t) + -^g= sinh( A /uJ " t) 
y(t) = y cosh( v /^ t) + -j$= smh(y/w^ t) 
z(t) — z + tw 

but in this case the condition (|14|) is satisfied if and only if xo = yo = and 
uq = vo = 0. This situation cannot occur. 

3. Review on the classical magnetic field on E 3 

As we have already said in Introduction, the best known example of magnetic fields in the 
Euclidean space E 3 is furnished by the 2-form F = dx A dy, corresponding to the Killing 
vector field £o = -§^- 

In this section we consider the Killing magnetic field = s Fq = s dx A dy, determined by 
the Killing vector field £ = s £o = s d z on E 3 , where s ^ is an arbitrary constant. We 
briefly describe its magnetic curves. 

The action of the Lorentz force on the vector space M 3 is given by: 

* f a = sd v , $^ = -sd x , = 0. 

Solving the Lorentz force equation 7" = $^(7^), we obtain the family of magnetic curves 
j s (t) = (x(t),y(t), z(t)), parametrized by 

' x(t) = Bi sin(st) + f cos{st) +x - f 
< y(t) = -Jcos{st) + J sin(st) +y + 2 f 
^z(t) = w t + z . 

Write the first Frenet equation 

7" = kN 

where k is the curvature and N is the normal of the curve. Using the equation (|4|) we obtain 
that the square of the curvature is 

K 2 = S 2 (l - Wl). 

Moreover, classical computations give the torsion r = sw$. 

Notice that even both the curvature k and the torsion r depend on the strength s, the ratio 
- does not. 

K 

We conclude with some comments: 

i) If Wo = the curvature is k = s and the torsion is t = 0, so the magnetic line is a 
(planar) circle. 

ii) If wo = ±1, then k = 0, r = ±s, so the magnetic curves are vertical lines. 

iii) In other cases the magnetic curves are circular helices. 
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4. Appendix 

In this section we present a Matlab program in order to compute, by numerical approxima- 
tion of the involved integrals, the parametrization of magnetic curve obtained in case I (i) 
from page [5] Since the curve is planar we consider <f> = 0. 



clear all 

7o7o7o Compute the integral 1(f) as a Riemann sum 

rho0=1.41; 
f _max=2 ; 
N=1000; 

L=(f_max-rho0~2)/N; 

for K=1:N+1 
a=0.001; 

b=rhoO~2+(K-l)*L; 

n=1000; 

h=(b-a)/n; 

k=0:n-l; 
x=a+k*h; 

f =1 . /sqrt (x . * (4-x . ~2) ) ; 

I(K)=h*sum(f ) ; 
J(K)=b; 

end 

•/.•/. \phi_0=0 

xx=sqrt ( J) ; 
y.yy=0*J; 

zz(l)=0; 
for K=1:N 

zz(K+l)=zz(K)-0.5*(I(K+l)-I(K))*J(K) ; 

end 

/o°/o the curve is planar 
plot(xx,zz, 'g-') 

text (0.25,-0. 75, '\rho_0=l. 41' , 'Color' , 'g') 
hold on 



Representation of the magnetic curves depending on the initial position: 
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P =0.001 \ 


j : 


p<f 1 \ 




P =1.2 \ 




P -1.3 




P -1.41 





1.9 1 ' 1 1 

0.5 1 1.5 



Figure 3. I (i) 

Using the previous Matlab program adapted to the example furnished at page [U and for 
the initial data ipo = and zq = 0, we can represent the corresponding magnetic curves: 

Magnetic curve in E 3 




Figure 4. x a = 1, y = 0, z = 0, u = 0, v = y 3 and w = ^-f^ 
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